CLAIMS 

What is claimed is: 

1 . A method of modeling faulting and fracturing in a subsurface volume of 
the earth comprising: 

(a) selecting a mode of definition of a subsurface model, said mode of 
definition selected from (i) an aerial mode wherein the model 
comprises a plurality of nodes in a horizontal plane interconnected 
to each other and to a substrate, (ii) a cross sectional mode wherein 
the model comprises a plurality of nodes in a vertical plane 
interconnected to each other and a substrate defining the edges of 
the model, and, (iii) a 3-D mode wherein the model comprises a 
plurality of nodes interconnected to each other and to a substrate 
defining the edges of the model; 

(b) defining said subsurface model including specifying material rock 
properties within the subsurface volume; 

(c) specifying an initial deformation pattern; and 

(d) using a dynamic range relaxation algorithm to find a force 
equilibrium solution for said subsurface model and said initial 
deformation pattern giving a resulting deformed model including 
fracturing. 

2. The method of claim 1, wherein selecting said mode of definition, defining 
a subsurface model, and specifying said initial deformation pattern further 



comprises using a graphical user interface. 
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3. The method of claim 1 , wherein said nodes are arranged in a grid that is 
one of (i) a triangular grid, and, (ii) a random grid. 

4. The method of claim 1 , wherein said nodes are interconnected by springs, 
and defining said subsurface model further comprises defining a normal force 
associated with each spring. 

5. The method of claim 4, wherein defining said subsurface model further 
comprises a substrate attachment force associated with each node that is attached 
to said substrate. 

6. The method of claim 1 , wherein said nodes are interconnected by beams, 
and defining said subsurface model further comprises defining a normal force and 
a shear force associated with each beam. 

7. The method of claim 1 , wherein specifying said initial deformation pattern 
further comprises performing a reconstruction based at least in part upon an 
observed large-scale deformation corresponding to said subsurface volume. 

8. The method of claim 7, wherein said reconstruction is a palinspastic 
reconstruction. 

9. The method of claim 7, wherein obtaining said initial deformation pattern 
further comprises: 

(i) obtaining a trial deformation pattern from said observed large scale 
deformations, 

(ii) applying an anticipate method to said model using said trial 
deformation, giving an approximate deformation result wherein said 
approximate deformation result is exclusive of fractures or faults; and 
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8 (iii) updating said trial deformation based on a comparison of said 

y approximate deformation result and said observed large scale deformation 

10 thereby giving said initial deformation pattern. 

1 10. The method of claim 7 further comprising conditioning said subsurface 

2 model thereby increasing the likelihood of said resulting deformed model 

3 including said observed large scale deformations, said conditioning including a 

4 weakening of bonds between adjacent ones of said plurality of nodes over at least 

5 a portion of the subsurface model. 

1 11. The method of claim 1 , wherein using the dynamic range relaxation 

2 algorithm further comprises applying said initial deformation model to said 

3 substrates in a plurality of steps, each step comprising a applying specified 

4 fraction of the initial deformation to said substrates and determining if any bonds 

5 between the nodes have been deformed beyond a breaking point and identifying a 

6 subset of the bonds that have been so deformed. 

1 12. The method of claim 1 1, wherein applying the dynamic range relaxation 

2 algorithm further comprises iteratively breaking the one bond of the subset of 

3 bonds that has been deformed the most and applying a relaxation algorithm to the 

4 remaining unbroken bonds. 

1 13. A method of simulating faulting and fracturing in a subsurface volume 

2 modeled by a plurality of interconnected nodes due to an initial deformation 

3 pattern applied to boundaries of said subsurface volume, the method comprising 

4 the following steps: 

5 (a) applying a fraction of said initial deformation to said boundaries; 

(> (b) using a dynamic range relaxation algorithm (DRRA) to find a force 
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equilibrium solution for said applied fractional initial deformation and 
identifying bonds between said plurality of interconnected nodes 
susceptible to breakage, wherein identifying bonds susceptible to breakage 
further comprises comparing a deformation of each bond in said force 
equilibrium solution to at least one predetermined breakage threshold 
associated with each of said bonds; 

(c) if in step (b) no bonds susceptible to breakage are identified, 
increasing said initial fractional deformation of said boundaries and 
increasing said fraction of said initial deformation and iteratively repeating 
steps (a) - (c ) until said fraction equals one; and 

(d) if in step (b), at least one bond is susceptible to breakage is 
identified, then breaking the one of the identified bonds susceptible to 
breakage whose deformation exceeds its at least one associated breakage 
threshold the most, and repeating step (b) 

14. The method of claim 13, wherein using a DRRA further comprises: 

(i) relaxing the plurality of interconnected nodes according to a single 
over-relaxation step to give a relaxed position of said plurality of nodes; 

(ii) identifying a first subset of the plurality of nodes that move further 
than a relaxation threshold; 

(iii) if said first subset of the plurality of nodes is empty, using said 
relaxed positions as said force equilibrium solution. 

15. The method of claim 14, wherein said first subset of nodes at step (iii) is 
not empty, and using the DRRA further comprises: 

I. relaxing sequentially each of the nodes in said first subset of nodes 
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4 and identifying a second subset of the plurality of nodes comprising those 

5 of the first subset of nodes, and each node connected thereto, that move 

6 more than a relaxation threshold ; 

7 II. after step II, interchanging the nodes in the first and second subset 

8 of nodes; and 

9 III. interatively repeating steps I and II until the first subset of nodes is 
1 () empty. 

1 16. The method of claim 13, wherein the connections between the plurality of 

2 interconnected nodes comprise a plurality of springs and at least one associated 

3 breaking threshold is an extensional breaking threshold. 

1 17. The method of claim 16, wherein the plurality of extensional breaking 

2 thresholds comprise one of (i) a Gaussian distribution characterized by a mean 

3 value and a standard deviation, and, (ii) a Weybull distribution. 

1 1 8. The method of claim 13, wherein the connections between the plurality of 

2 interconnected nodes comprise a plurality of beams and at least one associated 

3 breaking threshold comprises an extensional breaking threshold and a shear 

4 breaking threshold. 

1 1 9. The method of claim 1 8, wherein the plurality of extensional breaking 

2 thresholds comprise one of (i) a Gaussian distribution characterized by a mean 

3 value and a standard deviation, and (ii) a Weybull distribution, and the plurality of 

4 shear breaking thresholds comprise one of a Gaussian distribution and a Weybull 

5 distribution. 

1 20. The method of claim 13, wherein the subsurface volume further comprises 

2 a plurality of regions, each of said plurality of regions characterized by an 
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associated materials having material properties. 

21 . The method of claim 20, wherein said associated materials are selected 
from the group consisting of (i) salt, and, (ii) a rock. 

22. The method of claim 1 3, wherein the plurality of interconnected nodes 
constitute an aerial network and said connection between said interconnected 
nodes is selected from the group consisting of (i) springs, and, (ii) beams, and 
wherein the boundaries further comprise a plurality of substrate nodes, said 
plurality of substrate nodes attached to proximate nodes of the aerial network by 
the same type of connection as the connection between the interconnected nodes. 

23. The method of claim 13, wherein the plurality of interconnected nodes 
constitutes a 2-D cross section and wherein the boundaries comprise a plurality of 
discs, each of said plurality of discs experiencing at least one of (i) an attractive 
force towards, and, (ii) a repulsive force away from at least one of the plurality of 
interconnected nodes. 

24. The method of claim 13, wherein the plurality of interconnected nodes 
constitutes a 3-D network and wherein the boundaries comprise a plurality of 
spheres, each of said plurality of spheres experiencing at least one of (i) an 
attractive force towards, and, (ii) a repulsive force away from at least one of the 
plurality of interconnected nodes. 

25. The method of claim 23 further comprising checking a distance between 
pairs of said plurality of discs to a predetermined threshold after step (a) of claim 
1 3 and adding additional discs to said boundaries if said distance exceeds said 
predetermined threshold. 

26. The method of claim 24 further comprising checking a distance between 
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2 pairs of said plurality of spheres to a predetermined threshold after step (a) of 

3 claim 1 3 and adding additional spheres to said boundaries if said distance exceeds 

4 said predetermined threshold. 

1 27. The method of claim 1 8 further comprising applying a von Mises failure 

2 criterion. 

1 28. The method of claim 1 3 further comprising preconditioning the model to 

2 increase the likelihood of fracturing of the bonds between the plurality of 

3 interconnected nodes in the vicinity of specified locations. 

1 29. The method of claim 28, wherein said preconditioning further comprises 

2 reducing the predetermined breakage threshold of those of said bonds in the 

3 vicinity of said specified locations. 

1 30. The method of claim 29, wherein said specified locations are obtained by a 

2 geologic reconstruction based on observed large scale deformations. 

1 31. The method of claim 29, wherein said specified locations further comprise 

2 one of (i) piecewise linear curves, and (ii) piecewise linear surfaces. 

1 32. The method of claim 29, wherein reducing the strength of said bonds in the 

2 vicinity of the specified locations is based at least in part on a linear function of 

3 distance from said specified locations. 

1 33. A method of simulating deformation without faulting and fracturing due to 

2 an initial deformation pattern applied to boundaries of a subsurface volume 

3 modeled by a plurality of interconnected nodes, the method comprising: 

4 (a) defining a plurality of boundary nodes on a boundary of said 

5 subsurface volume wherein said initial deformation pattern is applied; 

<> (b) defining an initial and a final position for each of said plurality of 
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boundary nodes and a displacement there between; 

(c) determining a distance from each of the plurality of interconnected 
nodes to the final positions of the plurality of boundary nodes; 

(d) determining a displacement for each of the plurality of 
interconnected nodes as a combination of said displacement of said 
boundary nodes weighted by a weighting function related to said distance 
from each of the plurality of interconnected nodes and the final positions 
of the boundary nodes. 

34. The method of claim 33, wherein said weighting function includes an 
exponential factor related to said distance from each of the plurality of 
interconnected nodes and the final positions of the boundary nodes. 

35. The method of claim 34, wherein said weighting function further includes 
a matrix whose coefficients are obtained by solving a plurality of equations 
including said displacement of the boundary nodes wherein the plurality of 
equations is three limes the plurality of said boundary nodes. 

36. A graphical user interface (GUI) for displaying and manipulating a model 
of interconnected nodes for simulating fracturing and faulting in a subsurface 
volume of the earth, comprising: 

(a) a first module for presenting in a portion of a computer screen, a 
first graphical image representative of a plurality of interconnected nodes 
of the model; 

(b) a second module for defining material properties of the model 
defined by the plurality of interconnected nodes; 

(c) a third module for defining an initial deformation pattern applied to 



boundaries or the substrate of said plurality of interconnected nodes; 
(d) a fourth module for defining parameters of a simulation process 
including dynamic range relaxation algorithm for simulating a response of 
said model to said initial deformation pattern. 

37. The GUI of claim 36, wherein said first module is capable of presenting 
said first graphical image in one of (i) a planar view, and, (ii) a 3-D view. 

38. The GUI of claim 36 further comprising the input of a random number 
seed for the random number generator used for setting up at least one of (i) a 
geometry of said interconnected nodes in said model, and, (ii) breaking thresholds 
associated with links between pairs of interconnected nodes in said model. 

39. The GUI of claim 36, wherein said fourth module further comprises an 
editor for setting at least one of (i) a relaxation threshold for said simulation, (ii) 
an over-relaxation factor for said simulation, (iii) a maximum movement during 
said simulation, (iv) a time step for said simulation, (v) an angular relaxation 
factor lor said simulation, and, (vi) an angular over-relaxation factor for said 
simulation. 

40. The GUI of claim 36, wherein said third module further comprises an 
editor for defining a deformation that is at least one of (i) a localized extension 
button, (ii) uniform extension, (iii) uniform compression, (iv) a uniform right 
lateial shear, (v) a uniform left lateral shear, (vi) rotation, (vii) a deformation 
tegion in areal simulation mode or on the lowest plane in 3-D simulation, (viii) 
translation to a deformation region, and, (ix) a rotation to a deformation region. 

4 1 . The GUI of claim 36, wherein said second module further comprises a 
material editor for defining at least one region selected from the group consisting 
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of (i) a rock region in said model, and, (ii) a salt region in said model. 

42. The GUI of claim 4 1 , wherein the material editor further comprises 
defining, for each link associated with each pair of the plurality of interconnected 
nodes for the at least one region, properties selected from (A) an extensional 
breaking threshold for said link, (B) a shear breaking threshold for said link, (C) a 
linear force constant for said link, (D) a shear force constant for said link. 

43. 'I 'he GUI of claim 36, wherein said first module further comprises an 
editor for controlling the display of at least one of (i) faulting resulting from said 
simulation process, (ii) stresses resulting from said simulation process. 
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